Phase Separation Induced by Orbital Degrees of Freedom 
in Models for Manganites with Jahn- Teller Phonons 



00 

> 
o 



S. Yunoki, A. Moreo, and E. Dagotto 
National High Magnetic Field Lab and Department of Physics, Florida State University, Tallahassee, FL 32306 

(February 1, 2008) 

The two-orbital Kondo model with classical Jahn- Teller phonons is studied using Monte Carlo 
techniques. The observed phase diagram is rich, and includes a novel regime of phase separation 
induced by the orbital degrees of freedom. Experimental consequences of our results are discussed. 
In addition, the optical conductivity a{u!) of the model is presented. It is shown to have several 
similarities with experimental measurements for manganites. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 
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The properties of doped manganites are currently un- 
der much investigation due to the dramatic decrease in 
their resistivity when the spins order ferromagnetically 
by lowering the temperature T or applying a magnetic 
field jl|. This effect is caused by a metal-insulator tran- 
sition associated with the magnetic ordering. The ex- 
istence of the ferromagnetic (FM) phase is understood 
based on the double-exchange (DE) mechanism [|| . How- 
ever, experiments on manganites have revealed a compli- 
cated phase diagram that also includes charge-ordered 
and antiferromagnetic (AF) phases This rich struc- 
ture is beyond the DE ideas and a more refined approach 
is needed to understand these compounds. 

Since the 1950s the 1-orbital FM Kondo model for 
manganites has been widely studied. However, it is 
only recently that its computational analysis started, and 
surprises have already been observed In particu- 

lar, the transition from the undoped spin-AF regime to 
the spin-FM regime at finite hole-density occurs through 
phase separation (PS), instead of through a canted state 
as believed before. A growing body of experimental 
results indeed indicate the existence of PS in mangan- 
ites in agreement with the theoretical calculations. 
It is conceivable that PS tendencies above the FM criti- 
cal temperature Tf^^ could be important to explain the 
colossal magnetoresistance of doped manganites. 

However, in spite of its rich phase diagram the 1-orbital 
Kondo model is incomplete for a full description of Mn- 
oxides. For instance, dynamical Jahn- Teller (JT) distor- 
tions are also important Q, and a proper description of 
the recently observed orbital order [|l^ obviously needs 
at least two orbitals. Such a multi-orbital model with JT 
phonons is nontrivial [TT| , and the previous experience 
with the 1-orbital case suggests that a computational 
analysis is crucial to understand its properties. In ad- 
dition, it is conceptually interesting to analyze whether 
PS I^J^ exists also in multi-orbital models. 

It is precisely the purpose of this paper to report the 
first computational study of a 2-orbital model for man- 
ganites including JT phonons. The results show a rich 



phase diagram including a novel regime of PS induced 
by the orbital, rather than the spin, degrees of freedom 
(DOF). The Hamiltonian used here has three contribu- 
tions Hkjt = Hk + HjT + Haf- The first term is 



Hk = - ^ tab{c\a^C-.^ba- 
(ij)cra& 
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where (ij) denotes nearest-neighbor lattice sites, Jh > 
is the Hund coupling, a,b = 1,2 are the two e^-orbitals, 
the t2g-spins Si are assumed to be classical (with |Si| = 1) 
since their actual value in Mn-oxides (3/2) is large p2[ , 
and the rest of the notation is standard. None of the 
results described below depends crucially on the set {tab} 
selected . Throughout the paper the energy units are 
chosen such that tn = 1 in the i-direction. In addition, 
since Jh is large in the real manganites, here it will be 
fixed to 8 unless otherwise stated. Finally, the Cg-density 
(n) is adjusted using a chemical potential /i. 
The coupling with JT-phonons is through |^ 



iaba 



+ Q^'Kl (2) 



where Ql^ = -Qf - Qf\ and Qf = Qf = Qf\ 
These phonons are assumed to be classical. This approx- 
imation has been used and discussed in previous liter- 
ature where it was concluded that at temperatures 
of the order of the critical ones (room temperature) , or 
a sizable fraction of them, the use of classical phonons 
captures the important physics of the model ||lj]. Cer- 
tainly at very low-temperatures the quantum charac- 
ter of phonons is important, but this is not the range 
of temperatures explored in the present paper. Note 
that here T = 1/10 is about 200-300K g. Finally, a 
small coupling between the i2g-spins is needed to ac- 
count for the AF character of the real materials even 
when all La is replaced by Ca. This classical Heisen- 
berg term is Haf = ' Sj, where J' is fixed 

to 0.05 throughout the paper, a value compatible with 
experiments psl. To study Hkjt a Monte Carlo (MC) 
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technique was used. The trace over the Cg-electrons is 
carried out exactly using Hbrary subroutines for a fixed 
background of t2g-spins and phonons. This background 
is selected based on a Metropolis MC procedure [Q. 
The CPU time of the technique grows rapidly with the 
number of sites L, but the method has the important 
advantage that it does not have sign-problems at any 
T Finally, to analyze orbital correlations the pseu- 

dopin operator Ti = i Y^aab 4aa'^abCiba is used, while for 
spin correlations the operator is standard. The Fourier- 
transform of the pseudospin correlations is defined as 
Ti^l) = rEi,me''''^'~'"^Trn • Ti), with a similar defi- 
nition for the spin structure factor .^(q). 

Consider first the limit (n) = 1.0 (undoped mangan- 
ites). Fig. la shows T{q) and S{q) at representative mo- 
menta vs. A. For small A a large S{0) indicates a ten- 
dency to spin-FM order induced by DE (as in the quali- 
tatively similar 1-orbital problem at (n) = 0.5 The 
small values of T{q) imply that in this regime the orbitals 
remain disordered. When the couphng reaches Ad ^ 1.0, 
the rapid increase of ^^(Tr) now suggests that the ground 
state has a tendency to form a staggered (or "antiferro" ) 
orbital pattern, with the spins remaining FM aligned 
since 5(0) is large The existence of this phase was 
discussed before, but using multi-orbital Hubbard mod- 
els without phonons |jl^ . Our results show that it can be 
induced by JT phonons. As the coupling increases fur- 
ther, another transition at Ac2 ~ 2.0 occurs to a spin-AF 
orbital-FM state (5(7r) and r(0) are large). In this re- 
gion a 1-orbital approximation is suitable. Studying the 
spin and orbital correlations in real-space leads to the 
same conclusions as discussed here. 
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FIG. 1. (a) T{q) and S{q) vs A, working at (n) = 1.0, 
T = 1/75, Jh = 8, J' = 0.05, and in ID with 10 sites, {tab} 
correspond to set Ti (see [13]); (b) Same as (a) but for a 4^ 
cluster, T — 1/50, and hopping Ta {T4) in the y (x) direction. 
q — 0(7r) denotes (0,0) ((7r,7r)); (c) Same as (a) but for a 4'^ 
cluster, T — 1/50, the 3D hopping amplitudes of Ref. [13], 
and Jh ~ 00. q = 0(7r) denotes (0,0,0) ((tt, tt, tt)). 



limit where A and Jh are the largest scales, and using 
ti2 — ^21 = 0, til = t22 = t for simplicity. For parallel 
spins with orbitals split in a staggered (uniform) pattern, 
the energy per site at lowest order in t is — i^/A 0), 
where A is the orbital splitting. For antiparallel spins 
with uniform (staggered) orbital splitting, the energy is 

t^/2JH ( tV(2Jff + A)). Then, when A < 2Jh 

("intermediate" As), a spin-FM orbital- AF order domi- 
nates, while as A grows further a transition to a spin-AF 
orbital-FM ground state is expected. This reasoning is 
dimension independent, as the results for a 2D cluster in 
Fig. lb show. In 3D (Fig.lc) and Jh — 00 at least two 
of the regimes of Fig.la-b have been identified It 
was also observed that the behavior in Fig.la-b does not 
change when other sets {tab} are used [ll3|| . 

The next issue to be explored are the transport proper- 
ties in the three regimes at (n) = 1. The algorithm used 
here allows us to calculate real-time dynamical responses 
accurately, including the optical conductivity a{uj > 0), 
since all the eigenvectors in the fermionic sector for a 
given spin and phonon configuration are obtained ex- 
actly |^|. From the sum-rule, eg kinetic-energy, and the 
integral of (j{uj > 0), the ui = Q Drude- weight Dw can 
be obtained. In Fig. 2a, Dw is shown for several sizes. 
Dw vanishes at Ad signaling a metal-insulator transi- 
tion (MIT). Here the insulating phase is spin-FM and 
orbital- AF, while the metallic one is spin-FM and orbital- 
disordered |jl§]. The density of states (DOS) for A > Ad 
was also calculated and it presents a clear gap at the 
Fermi level. Although finite-size studies for D > 1 are 
difficult, the qualitative shape of Dw vs A on 4^ and 4^ 
clusters was found to be the same as in ID and, thus, it 
is likely that the MIT exists also in D = 2 and 3. 

Consider now the influence of hole doping on the 
(n) = 1.0 phase diagram. The first issue to be addressed 
is the stability of other densities as fi is varied. Fig. 2b 
shows (n) vs fi in the intermediate-A regime. It is re- 
markable that two regions of unstable densities exist be- 
low some critical temperature T^^ (roughly ^ 1/20, see 
inset Fig. 2b). Similar conclusions were reached using the 
Maxwell's construction Over 10^ MC sweeps at each 
H were needed for convergence near the unstable regions. 
These instabilities signal the existence of PS in the Hkjt 
model. At low-density there is separation between an (i) 
empty e^-electron band with AF-ordered i2g-spins and 
a (ii) metallic spin-FM orbital-FM phase. In the unsta- 
ble region near (n) = 1.0 PS is between the phase (ii) 
mentioned above, and (iii) the insulating spin-FM and 
orbital- AF phase described in Fig. la |20|. The driving 
force for this novel regime of PS are the orbital DOF, 
since the spins are uniformly ordered in both phases in- 
volved. Studying (n) vs /x, for A < Ad only PS at small 
densities is observed, while for A > Ac2 the PS close to 
(n) = 1 involves a spin-AF orbital-FA/I phase n . 



The three regimes of Fig. la can be understood in the 
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FIG. 2. (a) Dw vs A for several chains (T = 1/75); (b) 
(n) vs at A = 1.5, L = 10, and T = 1/40 (solid circles). The 
solid line is obtained from the Maxwell's construction. The 
triangles are results also at A = 1.5 and T — 1/40, but using 
14 sites and only 2 x 10* MC sweeps to show the appearance of 
hysteresis loops as in a first-order transition. The inset shows 
the T-dependence of the results at L = 10; (c) AE{n) (see 
text) vs (n) showing a negative curvature characteristic of PS. 
Results at large (open circles) and small (full squares) densi- 
ties are shown, on a L = 14 cluster with periodic boundary 
conditions and the couplings of (b). 

To confirm that the discontinuity in {n) vs /i corre- 
sponds to phase separation, in Fig. 2c the ground state 
energy is provided at several densities. The results indeed 
have the negative curvature characteristic of PS, both at 
large and small densities. To accommodate the two im- 
portant density regimes in the same plot, the energies in 
Fig. 2c are defined as AE{n) = E{n)—Eo{n), where E{n) 
is the actual ground state energy obtained as explained in 
Ref. 1^, and Eo{n) is a straight line (zero curvature) that 
joins the energies of the two (stable) extremal densities 
of both the low- and high-density regimes. Further con- 
firmation of the PS tendencies was obtained by inspect- 
ing the MC time-evolution of the density, as /i is varied. 
While at the critical /z frequent tunneling events were 
observed, at other /i's the time-evolution was smooth. 

Results in the limit Jh = oo using a L = 22 site cluster 
have also been obtained (see Fig. 3a). Once again, a dis- 
continuity in {n) vs /i was found at large and small den- 
sities, correlated with a negative curvature in the energy 
(not shown). Fig. 3a helps in clarifying that the plateaus 
at densities, e.g., between 0.3 and 0.5 in Fig. 2c are caused 
by the intrinsic discreteness of the clusters used here. As 
L grows, this fine structure disappears. See, e.g., the 
smoothly varying density between (n) ~ 0.25 and ^ 0.55 
in Fig. 3a. However, the discontinuities at large and small 
(n) remain strong, as they should if there is PS in the 
problem. Note that similar results as found in ID ap- 
pear also in studies of 2D systems (Fig. 3b). 




2.0 
1.5 



1 .0 1 1. 



;{b) 2D 












■ T=1/50 




0.5 ■ 








fij.; 




0.0 - 



_S-AF I 
O-F 1^ 
insulator i 



PS 

_S-F 
0-AF 
insulator,,, 

„ 



1 1 



AF 



■H 

S-F ^_ 
O-F 
metal 



PS 



S-F 
0-D 
metal 



w 



(c) 



-11.0 -10.0 -9.0 



-8.0 



1.0 O.i 



0.6 0.4 



0.2 0.0 



FIG. 3. (a) (n) vs fi in the limit Jh ~ oo, using A = 1.5, 
J' = 0.05, T = 1/30, and a 22-site chain; (b) (n) vs /i at 
A — 1.2 using a 4^ cluster and at T = 1/50. The two sets 
of points that produce the hysteresis loops are obtained by 
increasing and decreasing fj, using ~ 10* sweeps at each /i; (c) 
Phase diagram of Hkjt at Jh ~ 8.0, J' = 0.05, and using set 
Ti for {tab}- S-F and S - AF denote regimes with FM- and 
AF-spin order, respectively. O — D, O — F, and O — AF rep- 
resent states with disordered, uniform and staggered orbital 
order, respectively. PS means phase separation. 

Then, based on the information discussed thus far, sup- 
plemented by other MC measurements, the phase dia- 
gram of the ID Hkjt model is given in Fig. 3c. The 
metallic spin-FM region contains two regimes: one ferro- 
orbital ordered and the other orbitally disordered, as de- 
duced from the behavior of pseudospin correlations, the 
mean-value of the pseudospin operators, and the prob- 
ability of double occupancy of the same site with dif- 
ferent orbitals. The results are similar for several {tab} 
sets ||l^]. Our simulations suggest that the qualitative 
shape of Fig. 3c should be valid also in D = 2 and 3. 

Consider now (j(uj). Experimental studies for a variety 
of manganites reported a broad peak at w ~ leV (for 
hole doping x > 0.2 and T > Tf *^) ||l|,||. At room-T 
there is negligible weight near oj — 0, but as T is reduced 
the leT^-peak shifts to smaller energies, gradually trans- 
forming into a Drude response well below T^f*^. The 
finitc-cj peak can be identified even inside the FM phase. 
The coherent spectral weight is only a small fraction of 
the total. Other features at larger energies ^ 3eV involve 
transitions between the J/f -split bands and the 0-ions. 

In Fig. 4a, (j{uj) for the Hkjt model is shown at 
(n) = 0.7 and several temperatures near the unstable 
PS region of Fig. 3c (weight due to Jh split bands is not 
shown, but it appears at higher energy). Here the FM 
spin correlation length grows rapidly with the lattice size 
for T* < 0.05t, which can be considered as the "critical" 
temperature. Both at high- and intermediate-T a broad 
peak is observed at w 1, smoothly evolving to lower en- 
ergies as T decreases. The peak can be identified below 
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T* as in experiments |2l|,^. Eventually as the temper- 
ature is further reduced, (j{uj) is dominated by a Drude 
peak. The T-dependence shown in Fig. 4a is achieved at 
this A and {n) by a combination of a finite-cj phonon- 
induccd broad feature that loses weight, and a Drude re- 
sponse that grows as T decreases (for smaller As, the two 
peaks can be distinguished even at the lowest tempera- 
ture shown in Fig. 4a). The similarity with experiments 
suggests that real manganites may have couplings close 
to an unstable region in parameter space. In the inset, 
Dy^/^ vs T is shown. Note that Dy/ vanishes suggesting 
a MIT, probably due to polaron localization. Results for 
the 1-orbital case are smoother. 
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FIG. 4. (a) cr(Lj) at A = 1.0, (n) = 0.7, and L = 20, and 
several temperatures. The inset shows Dw vs T for both the 
Hkjt (circles) and the 1-orbital model (squares) of Ref. [4] 
(the latter at (n) = 0.65). Dw is normalized to its maxi- 
mum value at r = 0.01; (b) (^(ijj) vs uj parametric with (n) at 
A = 1.5, T = 1/10, and L = 16 (results for L = 10 are simi- 
lar). The inset shows the lower Jfr-split DOS at (n) = 0.93. 
In (a) and (b) a (5-function broadening e = 0.25 was used, as 
well as set Ti for the hopping amplitudes. 

A similar good agreement with experiments was ob- 
served working in the regime of the orbitally-induced PS 
but at a temperature above T^^ . Here the broad feature 
observed at high-T in Fig. 4a moves to higher energies 
(Fig. 4b) since A has increased. At the temperature of 
the plot the system is an insulator at (n) — 1, but as 
hole carriers are added a second peak at lower energies 
develop, in addition to a weak Drude peak (which car- 
ries, e.g., just 1% of the total weight at (n) — 0.61). This 
feature at high-T is reminiscent of recent experimental 
results where a two-peak structure was observed at 
room-T and several densities. Similar results were ob- 
tained on 4^ clusters. In Fig. 4b the peak at large-w is 
caused by phononic effects since its position was found 
to grow rapidly with A (A - 2A((Q(2)^ + Q(3)2)1/2^) 
corresponds to intersite transitions between Mn'^+ JT- 
split states. The lower energy structure is compatible 
with a Mn3+-Mn''+ transition pi. The inset of Fig.4b 



shows the DOS of the system. The two peaks above /j, 
are responsible for the features found in ajco). This in- 
terpretation is the same as given in Ref. |ll| at D = oo. 

Summarizing, here the first computational study of 
the 2-orbital Kondo model including Jahn- Teller phonons 
was reported. The phase diagram includes regions of 
phase separation both at large and small e^-density. One 
of them corresponds to a novel regime where PS is driven 
by the orbital, rather than spin, degrees of freedom. 
Coulomb interactions will break the large regions in- 
volved in PS for the pure Hkjt model into small islands 
of one phase embedded into the other. X-ray diffraction 
measurements should observe a coexistence of patterns 
characteristics of the two extreme stable phases, if the 
PS is static. If the process is more dynamical the Bragg 
peaks should be broad. Recalling that our spin-FM phase 
at (n) — 1 and D — 1 , 2 is also compatible with A- 
type spin-AF order, several experimental results |^ are 
in agreement with the tendencies to PS discussed in this 
paper. In the limit of small (n), the PS observed here 
could be transformed by long-range Coulomb interactions 
into a charge-ordered state. Note that the PS observed 
in electron-doped Sr2Mn04 ||^ is also compatible with 
our results Q. 
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